Nuclear level statistics: extending the shell model theory to higher temperatures 
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The Shell Model Monte Carlo (SMMC) approach has been applied to calculate level densities and 
partition functions to temperatures up to ~ 1.5 — 2 MeV, with the maximal temperature limited 
by the size of the configuration space. Here we develop an extension of the theory that can be used 
to higher temperatures, taking into account the large configuration space that is needed. We first 
examine the configuration space limitation using an independent-particle model that includes both 
bound states and the continuum. The larger configuration space is then combined with the SMMC 
under the assumption that the effects on the partition function are factorizable. The method is 
demonstrated for nuclei in the iron region, extending the calculated partition functions and level 
densities up to T ~ 4 MeV. We find that the back-shifted Bethe formula has a much larger range of 
validity than was suspected from previous theory. The present theory also shows more clearly the 
effects of the pairing phase transition on the heat capacity. 

PACS numbers: 21.10.Ma, 21.60.Cs, 21.60.Ka, 05.30.-d 



I. INTRODUCTION 



Consistent calculations of nuclear partition functions and level densities are required in nuclear astrophysics. The 



■ nuclear partition functions at a finite temperature are needed to understand pre-supernova collapse |1( and for the 
calculation of stellar thermal reaction rates Q ■ Similarly, level densities determine the statistical neutron and proton 

. capture cross-sections in nucleosynthesis Q. 

The Fermi gas formula (also known as the Bethe formula) Q| for the nuclear level density can be derived from the 
, partition function of non-interacting nucleons [5j. Correlation effects, e.g., shell, pairing and deformation effects, are 
usually accounted for through empirical modifications of this formula. In particular, the back-shifted Bethe formula 
0, Q is found to describe well the experimental data if its parameters, the single-particle level density parameter a 
and back-shift parameter A, are fitted for each nucleus H,I3- However, these phenomenological approaches cannot be 
reliably used in nuclei for which experimental data are not available, so there is a need for microscopic theories that 
include interaction effects. In mean field theory, pairing can be treated in the BCS approximation, as was done in 
Ref. [Toj . However, shape fluctuations are likely to be significant and they require a theory beyond the mean field, e.g., 
the static path approximation |ll| . A theory that includes quadrupole fluctuations in the static path approximation 
was proposed in Ref. [l2T |. Like the mean field theory, the static path approximation allows calculations in large 
spaces. This was exploited by the authors of Ref. 0] who presented results for large single-particle spaces including 

■ an approximate treatment of unbound states. However, as implemented there, the theory includes neither pairing 
correlations, which are certainly of similar importance, nor multipoles of the interaction beyond quadrupole. 

A more systematic approach to calculate correlation effects is to use the interacting shell model. In this approach, 
we define a many-particle configuration space and treat in full the effective interaction within that space. Here the 
shell model Monte Carlo (SMMC) method has rather favorable computational properties, scaling with the number 
N of single-particle orbitals as iV 4 , compared to direct diagonalization methods, which scale exponentially. 1 Even 
so, most applications of the SMMC have so far been limited to a single major shell. For medium mass nuclei, this 
restricts the reliability of the calculated partition functions (and associated level densities) to temperatures below 
~ 1.5 -2 MeV Hi EI Ei- 

In this work, we will first explore the model space truncations in an independent-particle approximation, extending 
the single-particle space to include both the bound states and the continuum. This was studied within a semiclassical 
approximation in Ref. poj . We will develop the fully quantum- mechanical theory that includes the continuum 2 in 



1 Another m ethod in the literature is the use of moments of the Hamiltonian to expand the level density in a sum of partitioned Gaussians 
11 1 Il4l (Li. For a similar approach using binomials instead of Gaussians see Ref. 1 1(1 . 

2 In Ref. Il2l the continuum was treated approximately by enclosing the nucleus in a sphere whose radius R is slightly larger than the 
size of the nucleus. Our treatment is exact and corresponds to the limit R — * oo. 
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Section^ We find that in tightly bound nuclei, the continuum contribution is negligible even at high temperatures 
(e.g., T ~ 4 MeV for 56 Fe), while in weakly bound nuclei the continuum contribution can become significant already 
at low temperatures. In either case the continuum contribution can be well approximated by considering just the 
contribution of the narrow resonances, treated as broadened bound states. 

The remaining task is to include interaction effects, as discussed in Section lTTTl It is necessary to distinguish between 
the long-range part of the interaction, acting in the valence orbitals, and the short-range part, involving highly excited 
configurations. The short-range part is beyond the scope of the theory, and must be taken into account in defining 
an effective long-range interaction, but the latter can be accurately treated within the SMMC approach. At higher 
temperatures, long-range interaction effects are weak and corrections to the entropy from the larger model space 
can be included in a single-particle approach. We shall see that at high temperatures, the interaction effects on the 
partition function scale as 1/T, and the two temperature regimes can be joined smoothly. We find that the extended 
partition function and its associated level density are well described by the back-shifted Bethc formula up to T ~ 4 
MeV. This value represents an upper limit for our underlying assumption that the single-particle potential is fixed 
and independent of temperature. 

The basic object we study is the canonical partition function defined by 

Z((3) = Tre-P" , (1) 

where f3 is the inverse temperature, [3 = and the trace is taken at a fixed particle number (i.e., fixed number 
of protons and neutrons). We shall find it convenient to measure the energy with respect to the ground energy Eq, 
defining an excitation partition function Z' as 

Z' = Ze 0Ea . (2) 

Thus at zero temperature, the excitation partition function is equal to the degeneracy of the ground state N gs , 
Z'(T = 0) = N gs . 

The conversion from the partition function to the level density is discussed in Section IIVI While the results of the 
theory can be given as numerical tables of Z(T) and the corresponding level density p(E), it is also useful to express 
the results in terms of the parameters of the back-shifted Bethe formula, which provides a good description at not too 
low temperatures. The inclusion of a back shift requires a particular parameterization of Z', as discussed in Section 

eh 



II. ORBITAL TRUNCATION EFFECTS IN THE INDEPENDENT-PARTICLE APPROXIMATION 

When the temperature is larger than about a third of the nucleon separation energy, the unbound nucleon states 
can no longer be ignored. To develop a theory that includes the continuum, we first need to define precisely how we 
separate the contributions of the free nucleons and the nucleus to the partition function. If the system is enclosed 
in a box, the free nucleon partition function is proportional to the volume of this box. The contribution of the 
continuum to the nuclear partition function is found by the explicit subtraction of the free-particle contribution, and 
is independent of the volume of the box, if it is large compared to the size of the nucleus. The free nucleon partition 
function we use here is the one obtained when the nuclei are ignored entirely. 3 

In the independent-particle approximation, the single-particle spectrum is calculated in an external potential well, 
representing the nuclear mean-field potential. Assuming the well to be spherically symmetric, the bound state energies 
e n ij depend only on the principal quantum number n, orbital angular momentum I, and total spin j. For the positive 
energy continuum, we need the scattering phase shifts <%(e) where e is the (positive) energy of the continuum states. 
The change Sp in the single-particle continuum level density in the presence of the potential V is found by subtracting 
the free particle level density po from the total density p: 

Sp(e) = p{e) p (e) = £(2j + 1)^ . (3) 

7r de 

h 

We denote the many-body grand canonical partition function in the independent-particle approximation by Z^ c , 



3 We also ignore the Coulomb effects on the partition function, which are small but which immediately introduce the complications of 
plasma theory. 
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where the subscript stands for "single-particle." It is given by 

/•OC 

\nZg c (j3, n) = V(2j + 1) ln[l + e H»(«««-iO] + / de6p(e) ln[l + e^ 6 ^] , 



(4) 



nlj 



where [i is the chemical potential. 4 A more convenient form, avoiding the phase-shift derivative, may be obtained by 
integrating the second term in the above equation by parts. The contribution from the integral endpoint at e = is 
evaluated using Levinson's theorem, 5y (e = 0) = nijir, where nij is the number of bound Ij states. The result is 



jr />oo 

de-± ln[l + e-^-^} = -n W 7rln(l + e^) +/3 / de %(e)/(e) , 
o " e Jo 



(5) 



where /(e) = [1 + e^ e 1 is the Fermi-Dirac occupation. The formula for the partition function can then be 
rewritten 



lnZ s G p c (/3, M ) = ^(2 J + l) £ln 



1 + e 



-P(e nlj -ii) 



l + e^ 



& ««(e)/(e) 



(6) 



This form is easier to use computationally because it avoids numerical derivatives. 

To calculate the partition function of a specific nucleus, we need to transform the grand canonical partition to 
the canonical partition function at fixed proton and neutron numbers. In the saddle point approximation, the 
corresponding correction is estimated from the particle-number fluctuations in the grand-canonical ensemble, giving 
for the canonical Z(f3) 



Z((3) « (27r((AA0 2 >r 1/2 Z GC e-W 



(7) 



or 



bxZ' w \nZ GG + (3E - (3(iN - - In (2tt((AA^) 2 )) 



(8) 



Here ((AiV) 2 ) is the variance of the neutron or proton number and Z' = Ze^ E °. In Eq. (JSJ, fi is calculated from 
(N) = d In Z GC I da — N (where (N) is the average particle number in the grand-canonical ensemble and a = /3fx), 
and the variance is found from ((AiV) 2 ) = d 2 In Z GC /da 2 . In the independent-particle approximation these quantities 
are 



dlnZ GG „ 

v = - =E( 2 j +1 ) 



da 



(9) 



, d 2 lnZ GG _ 

(f^ a > = -^=E( j y+ 1 ) 



E/ny(l-/ny) + 1 f°^/(l 
„ ^ Jo de 



f) 



To avoid numerical derivatives, we can again use integration by parts and rewrite Eqs. © and (|1U[) as 



JV = E(2J + 1) 



and 



((atv) 2 ) = E(2j + 1) { E t/^-C 1 - - /(oJf 1 - /(°))1 

Ij V n 



+£ / de%(e)/(l-/)(l-2/) 

7T 



(10) 



(11) 



(12) 



4 In practice, Eq. J1J includes two separate contributions from protons and neutrons with different chemical potentials fi p and fj,„ 
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FIG. 1: The logarithm of the canonical excitation partition function Z' sp (where the energy is measured with respect to the 
ground-state energy) calculated for 56 Fe is shown for the following single-particle spaces: all bound states plus continuum as 
described by Eqs. 0, JHJ, an d H'2t (solid line); bound states and the narrow resonances listed in Table I (dashed-dotted line); 
and bound states only (dotted line). The corresponding calculation with the space truncated to the fpg major shell is shown 
as the dashed line. The intercept at T = is IniVgs, where N gB = 168 is the number of states in the {fy 2 )~ 2 {P3/2) 2 shell 
configuration. 

To illustrate the method, we apply it to nuclei in the iron region. For the mean-field potential we use a Woods-Saxon 
well V(r) plus a spin-orbit interaction 

^) + A la ^s)r^, (13) 

where V(r) = Vq [1 + exp(r — Rq)/o\~ X with Rq = tqA 1 / 3 . We choose the parameterization of Ref. where Vq — 
—51 + 33 N ^ z MeV, a = 0.67 fm, tq — 1.32 fm, and Xe s — —0.44. We have calculated the single-particle spectrum 
and scattering phase shifts in this potential well, and used them to compute the canonical partition function from 
Eqs. © an d (|12|) . The solid line in Fig. ^ shows the logarithm of the canonical excitation partition function 
lnZg p for the nucleus 56 Fe when the full single-particle space (bound states plus continuum) is taken into account. 

TABLE I: Neutron orbital energies and resonance energies (in MeV) in 56 Fe, calculated with the Woods-Saxon potential of 
Ref. 0. 



Orbital type 




e ( 5B Fe) 


e ( 66 Cr) 


bound 


/r/2 


-12.68 


-10.2 




P3/2 


-8.65 


-6.7 




Pl/2 


-6.56 


-5.1 




h/2 


-6.34 


-5.2 




99/2 


-3.13 


-1.7 




db/2 


-0.41 






Sl/2 


-0.26 




resonance 


^5/2 




0.35 




97/2 


4.95 


4.67 




hu/2 


6.22 


6.62 



We next investigate various truncations within the independent-particle model. Of particular interest is the trun- 
cation to the same model space that is used in the SMMC approach. Most of the SMMC calculations so far have been 
restricted to the bound states in a single major shell. For example calculations in the iron region were restricted to 
the complete fpgg/2 major shell, which includes the active orbitals /7/2,P3/2>Pi/2j /s/2 and 39/2 ■ The orbital energies 
[determined in a Woods-Saxon well with spin-orbit interaction Ijl3(l ] are given in Table I. 
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FIG. 2: Similar to Fig. □ but for the nucleus 66 Cr. 



The partition function calculated in this truncated space (in the independent-particle approximation) is compared 
with the full space calculation [Eqs. ||BJ|, iJSJ and l|12|) ] in Fig. ^ One sees that the truncation to a single major shell 
becomes problematic at temperatures above ~ 1.5 MeV. 

To assess the importance of the continuum, we also show in Fig. ^ the result of the calculation keeping all bound 
states but neglecting the continuum integral. One sees that the continuum contribution is not very significant for a 
tightly bound nucleus such as 56 Fe. 

However this situation changes in a nucleus with a small neutron separation energy. In Fig. [2] we compare the 
calculations in the full single-particle space and in the truncated space that includes all bound states for the nucleus 
66 Cr. Here the separation energy of the highest occupied orbital is 1.7 MeV. One sees that the continuum effects 
become significant for T > 1 MeV. 

Rather than ignoring the continuum completely, one could also consider the approximation of including only narrow 
resonances, treating them as broadened discrete states. Indeed, the derivative of the scattering phase shift <5y in the 
vicinity of a narrow resonance follows a Lorcntzian 



de 



( e -^.) 2 + (iy2) 2 



(14) 



where ey and Fy- are the centroid and width of a resonance with quantum numbers Ij. If the resonance width Fy is 
much smaller than T, the integral in |0J can be evaluated to give a contribution to In Zf^ as though the resonance 
were a discrete state with energy ey- . 

We identify the resonances from the energy dependence of the phase shifts <$y in different channels Ij, and take their 
energies where the phase shifts ascend through ir/2. The channels with resonances and their energies are displayed 
in Table I. The partition function calculated using bound states and narrow resonances is shown as the dotted line in 
Fig. |2 We see that the resonance approximation is good at least to T ~ 3.5 MeV. 

We now summarize the validity of the various approximations, requiring that the theory be accurate up to tem- 
peratures of T ~ 4 MeV. For nuclei along the stability line, it is only necessary to include bound states. If one goes 
to nuclei that have much lower separation energies, the continuum contribution can be significant but the resonance 
approximation is likely to be adequate. It should be emphasized throughout this section we have assumed a fixed 
mean-field potential. Above T ~ 4 MeV, that approximation becomes unreliable. 



III. INTERACTION EFFECTS 



We next discuss interaction effects on the partition function. The ground-state binding energy is of course very 
sensitive to the interaction, but this will not be immediately visible in the excitation partition function. Let us denote 
by Z' v tr the partition function calculated with interactions but in a truncated single-particle space. We argue that the 
correction due to a larger model space will be largely additive in the logarithm of the partition function, so that the way 
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to combine a small-space calculation of interaction effects with a large-space calculation of the independent-particle 
partition Z' sp is with the formula 

In Z' v = In Z' v tr + In Z' sp - In Z' sptr , (15) 

where Z' sp tr is the excitation partition function (in the absence of interactions) in the same truncated single-particle 
space in which Z' v tr is calculated. 

Eq. i|15|) cannot be justified rigorously, but we can motivate it in several ways. It is certainly true that the interaction 
corrections approach zero at high temperature. From finite-temperature many-body perturbation theory, it is seen 
that the leading corrections to the independent-particle Hamiltonian H are additive in the logarithm of the partition 
function, as is assumed in Eq. Ijl5|l . The explicit formula is 

rP i rP rP 

lnTr cxp(-/3X)-lnTr exp(-(3K ) = - / rfr<Vr(r)) + - / dr x dT 2 (T T V I (r 1 )V I {T 2 )) conn + . . . , (16) 

Jo z Jo Jo 

where K = H — /iiV, Kq = Hq — fiN, Vi is the two-body interaction in the interaction picture (with respect to Kq), T T 
denotes time ordering and (. . .) denotes a thermal average with respect to Kq. Each integration in (|16fl gives a factor 
1/T. In cases where analytic expressions can be derived (e.g., hard core gas and Coulomb gas), it may be verified 
that the interaction correction falls off as T" 1 or more strongly at high T. At low temperatures, the extended space 
is not thermally occupied and only plays an indirect role, renormalizing the interactions in the smaller space. Since 
the SMMC is applied with renormalized interactions that are appropriate for the smaller space, there is no additional 
correction when the larger space is considered explicitly. With the correct limiting behavior at both low and high 
temperatures, we have a better theory for the complete temperature dependence. 

We now give a brief summary of the SMMC calculations to be used for Z' v tr in Eq. (|15f) . The interaction is taken to 
be the same as in Ref. |17| . It is separable and surface-peaked, acting between orbitals of a major shell (here the fpgg/2 
orbitals) . The strength of the surface-peaked interaction is determined self-consistently and renormalized appropriately 
for each multipole. In addition, there is a monopole pairing interaction, whose strength is determined from odd-even 
mass differences. The single-particle Hamiltonian is the same as in Section^ The outputs of the SMMC computation 
are canonical expectation values of various observables. For our purposes here, the most important quantity is the 
canonical thermal energy E((3), calculated as the expectation of the Hamiltonian with number-projected SMMC 
densities. The partition function is then found by integration 

rP 

InZ(p) = lnZ(O) - / df3'E{(3') . (17) 
Jo 

Here Z(0) is the partition function at T — > oo and is equal to the total number of many-particle states in the model 
space. In addition, we need to determine the ground-state energy Eq in order to find the excitation partition function, 
Eq. J2J). This is done by extrapolating the calculated E((3) to large (3. 

The SMMC results for the nucleus 56 Fe are shown in Fig. [21 with \yiZ' v tr plotted as a function of T (squares). The 
curve starts out flat, changing to a linear increase over some range of T, and reaching a plateau beyond the highest 
temperature plotted. The initial flat behavior is associated with the gap in the energy spectrum between the ground 
state and the first excited state, while the linear regime corresponds to the Bethe formula. The saturation is due to 
the truncation of the space, and its onset marks the limit of validity of the truncated SMMC calculation. We next 
repeat the SMMC calculation with the two-body interaction turned off. The results are shown by the circles in Fig. [3] 
Notice that the difference between the two partition functions becomes small as the temperature increases, confirming 
the analysis in an earlier part of this section [see in particular Eq. (|1GII ]. 

Before calculating the partition function in the extended model space, we also compare in Fig. two ways of 
calculating the non-interacting partition function, from Eq. (JSJ (dashed line) or from the SMMC with the interaction 
turned off (squares). In the latter method, particle- number projection is carried out exactly, while in Eq. (jHJ this is 
done in the saddle point approximation. One can see that the differences are insignificant for our purposes. 

We now combine the various terms in Eq. I|15fl to get our extended-range partition function. The result is shown in 
Fig. 0]by the solid squares. For comparison, the SMMC result (open squares) and the full space independent-particle 
model result (dashed line) are shown as well. Remarkably the extended In Z' is seen to be close to a linear function 
of T over most of the range, unlike either of its constituents. In the saddle point approximation, the dominating term 
in In Z' is linear in T with whose slope is the parameter a in the level density formula. However, to compare with the 
back-shifted Bethe formula [see Eq. <|27|l in Section HVj . more care is needed in treating the sublinear terms. Starting 
from the Darwin-Fowler expression for the partition function of non-interacting fermions [see Eq. (2B-9) in Ref. @], 
and including the proton and neutron number fluctuation terms from Eq. © , one can derive the following expression 
for Z 1 , 



In Z' «aT- ln(6aT/7r) , 



(18) 



N 10 




FIG. 3: The excitation partition function of 56 Fe using the SMMC method in the truncated space (fpgg/2)- The partition 
function Z' v tr in the presence of the residual interaction (squares) is compared to the SMMC partition function Z' sPitr in the 
absence of interactions (circles). The dashed line is the partition function of the independent-particle model calculated in the 
saddle point approximation using Eqs. © and il l 21 . 




T(MeV) 



FIG. 4: Extended partition function for 56 Fe, calculated from Eq. 1151 shown as solid squares. The solid line is a fit to 
Eq. (O with a = 5.87 MeV -1 and A = 1.03 MeV. Shown for comparison are the SMMC (open squares) and the full space 
single-particle (dashed line) results. 



where a — Tt 2 g/3 and g is either the proton or neutron single-particle level density at the Fermi energy (assumed to 
be equal). Within the saddle-point approximation, Eq. 118(1 is equivalent to the simple Bethe formula for the level 
density (see Section llVjl . As mentioned in the introduction, a better phenomenological description of the measured 
level density is often obtained by shifting the ground state energy by an amount A in the Bethe formula. We show 
in Appendix A that the corresponding modification of the partition function is to introduce a third term in Eq. 118(1 . 

In Z' ^aT - ln(6aT/7r) - A/T . (19) 

The solid line in Fig.0]is a fit of 1(19(1 to In Z' for 1 < T < 4 MeV with a = 5.87 MeV -1 (corresponding to the expression 
a = A/K with K = 9.5 MeV) and A = 1.03 MeV. The fit does remarkably well in describing the T-dependence of 
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T(MeV) 

FIG. 5: The extended entropy of 56 Fe (solid squares). The solid line is a fit to Eq. witn a = 5.82 MeV -1 . The open 
squares describe the SMMC entropy. 

\nZ' for temperatures between T ~ 1 MeV and T ~ 4 MeV (the rms of the deviation from the calculated hxZ 1 is 
0.05). The functional form i|19|) of InZ' is the essence of the (back-shifted) Bethe level density formula; the result 
observed here implies that this formula is useful to considerably higher temperatures than was apparent in earlier 
calculations. At high temperatures, the dominant term in ljl9|l is linear in T. However, the other two terms are 
necessary to obtain a value of a that is comparable to the value extracted directly from the level density, which will 
be the subject of Section IIVI 

In the independent-particle model, InZ' is approximately linear in T only for T > 2 MeV with a slope of a ~ 5.25 
MeV -1 , i.e., K ss 10.7. This value of K is smaller than the Fermi gas value 5] K 15 but still larger than the value 
we found above when correlations are taken into account. 

The canonical entropy can be calculated from S — —dF/dT where F = —T\n.Z' is the canonical free energy. Using 
the empirical formula (|19J) . this entropy is given by 

S w 2aT - ln(6aT/7r) - 1 , (20) 

and is independent of the back-shift parameter A. Thus the parameter a can be determined by a single parameter 
fit of the extended entropy to Ij20|l . The leading order term in the entropy is linear in T and the last two terms in 
(12011 originate in the particle number fluctuations. The microscopic calculation of the extended entropy can be done 
without taking numerical derivate with respect to T, as is explained in Section ll VI Fig. [S] shows a fit of the empirical 
formula i|20|) (solid line) to the extended entropy in the range T ~ 1 — 4 MeV (symbols) with a = 5.82 MeV -1 (the 
rms deviation is rs 0.3). This value of a is in agreement with the value found from In Z' . The fit is good and only at 
T ~ 4 MeV, do we start to observe a small deviation. 

We have demonstrated our method for an even-even nucleus 56 Fe, but it should be equally applicable to odd-even 
and odd-odd nuclei. As an example we show in Fig. the logarithm of the excitation partition function for 57 Fe. The 
solid line is a fit to {Eft with a = 6.04 MeV" 1 and A = -1.29 MeV in the range 1 < T < 4 MeV (the rms deviation 
is w 0.08). A qualitative difference is that unlike the case of 56 Fe, \nZ' does not approach zero at low temperatures. 
This is due to the spin degeneracy of the ground state of an even-odd nucleus and the much smaller energy gap to 
the first excited state. 



IV. LEVEL DENSITIES 

We now determine the level densities associated with the extended partition function of the previous section. The 
level densities are calculated from the partition function as in Refs. IT] and [3 using the saddle-point formula 



p{E) = {2Trf3- 2 C)- 1/2 e s , 



(21) 
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FIG. 6: The extended partition function (solid squares) for the odd-even nucleus 57 Fe is compared with the SMMC partition 
function (open squares). The solid line is a fit to the back-shift formula (1191 with a = 6.04 MeV -1 and A = —1.29 MeV. 



where S = In Z+/3E = In Z' +fl(E — Eq) is the canonical entropy and C = —f3 2 dE/df3 is the canonical heat capacity. 
The thermal energy in (|21|l is a function of (3 determined by 



E-Eq 



9 In Z' 

d(3 



Using the expression (|15|) for the extended partition function and Eq. (|22|l . we have 



E — En 



E 



[E sp — (E sPttr + E t 



) 

core / 



(22) 



(23) 



where E% n 



E° 

sp 



E% ttr is theT: 



energy of the core (i.e., the shells below the valence shell). The notation used 



here for the various energies follows the same notation used for the partition functions in (|15f) and the subscript or 
superscript "0" denotes the ground state in the corresponding model space. 
Similarly, we find for the extended entropy and heat capacity 



S — S Vj tr 
C = C v tr 



(24) 



Eq. (|23J) is used to determine the relation between the excitation energy and j3 in the calculation of the extended 
level density. When the SMMC approach is used to calculate interaction effects in the truncated space, E Vytr is 
calculated directly from the expectation value of the Hamiltonian and E® tr is found by extrapolation to (3 — > oo. The 
partition Z v ^ r (and therefore the entropy S v j r ) is found by integration using (|17|) . while the heat capacity C Vf tr is 
found by a numerical derivative of the energy. 

It remains to determine the respective quantities in the independent-particle model (both in the full and truncated 
spaces). Here there is a technical complication in that the needed quantities are canonical. For example, we need to 
calculate the canonical energies E sp and E sp ^ r - In principle, one can use l|22|l where Z is calculated in the independent- 
particle model. However, since Z is a canonical partition function, the partial derivative with respect to must be 
evaluated at constant particle number N rather than at constant a (a = 0fi) . The explicit expressions for E sp and C sp 
in terms of logarithmic derivatives of the grand-canonical partition function are derived in Appendix B. In Appendix 
C we express these logarithmic derivatives in terms of the corresponding single-particle spectrum and scattering phase 
shifts. 

The extended level density of 56 Fe is shown in Fig. [7| (solid squares) as a function of excitation energy, E x = E — Eq. 
For comparison, the SMMC level density calculated in the truncated space (i.e., fpgs/2 shell) is shown by open squares, 
while the level density calculated in the independent-particle model (no truncation) is shown by the dashed line. 

To derive the back-shifted Bethe formula for the level density, we evaluate the inverse Laplace transform of Z' in 
Eq. H18|) in the saddle point approximation (separating out the particle number fluctuation term as a pre-exponential 
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FIG. 7: Level density of 56 Fe, calculated with the extended partition function (solid squares). The solid line is a fit to the 
back-shifted Bethe formula I27H and describes well the extended level density in the full excitation energy range shown in the 
figure. The SMMC level density in the truncated fpgg/2 shell is shown by open squares and the independent-particle model 
level density is shown by the dashed line. 



factor). We find 



H V ' V dP) V6aT/ 12 
where the temperature is determined by the excitation energy E x = E — Eq through 

E x w aT 2 + A . 



-1/2 



(25) 



(26) 



Relation 1)26(1 generalizes the usual Bethe relation E x = aT 2 to include an empirical offset A to the ground-state 
energy. This offset originates in correlation effects. 

When expressed in terms of the excitation energy, Eq. 1251) is just the back-shifted Bethe formula 



P(E X ) 



-1/4 



12 



(E x - A)~ 5 / 4 exp[2\/a(£ , x - A)] 



(27) 



This formula has two fit parameters, the single-particle level density parameter a and the back-shift parameter A. 
Fitting H27fl to the calculated extended level density in the range 5 < E x < 60 MeV, we find a = 5.9 MeV -1 and 
A = 1.35 MeV (with a x 2 P er degree of freedom below 1), in overall agreement with the values found in Section ITTT1 
by fitting ((19(1 to \nZ'. These values of a and A are also in close agreement with the earlier values found by fitting 
to the SMMC level density in the energy range 5 < E x < 16 MeV, a = 5.85 MeV" 1 and A = 1.37 MeV. However, 
now the back-shifted Bethe formula 1(27(1 describes the extended level density up to much higher values of E x (~ 90 
MeV). Thus the validity of this formula is now extended to significantly higher excitation energies than previously 
thought. Furthermore, the specific values of the parameters a and A found in previous SMMC calculations [18J are 
approximately valid in this extended energy range. 

Fig. |S| shows the level density of the odd- A nucleus 57 Fe, where we find in the extended range a = 6.05 MeV" 1 and 
A = —1.3 MeV. These values are comparable to the values determined from InZ' in Section ITTll 

The good fits we find to the back-shifted formulas ((19(1 and ((27(1 , for the logarithm of the partition function and the 
level density, respectively, imply that the excitation energy can be approximated by the quadratic expression 1(26(1 for 
not too low temperatures. 



5 A temperature of T=4 MeV for 56 Fe corresponds to E x & 86 MeV. 
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FIG. 8: Extended level density (solid squares) in comparison with the SMMC level density (open squares) for 57 Fe. The solid 
line is a fit to the back-shifted Bethe formula 1271 , 

Experimentally, the statistical properties of nuclei at higher temperatures have been deduced from heavy ion reac- 
tions. If the nuclei are equilibrated when they decay, there is a direct relation between the kinetic energy distributions 
of the decay products and the temperature of the daughter nucleus. Experimental results have been reported indicat- 
ing that the quadratic dependence of E x on temperature needed to be modified at high temperature |21| . However 
in other experiments only a slight deviation was found |22| . These latter experiments concluded a weak decrease of a 
with temperature and even a constant a could not be excluded, in agreement with our results. 



The back-shifted parameterization works well for temperatures T 1 MeV, but not at lower temperatures. The 
deviations can be clearly seen in the canonical heat capacity as a function of the temperature. There can be a peak 
in this function that is associated with the pairing phase transition, but the back-shifted parameterization permits 
only a monotonic function rising linearly with temperature. Recently, experiments in rare-earth nuclei have been 
reported in which effects of pair breaking can be seen in the canonical heat capacity |2j 0] ■ The authors measured 
the level density, constructed a canonical partition function from their data, and used it to deduce the heat capacity 
as a function of temperature. 

Such low temperature characteristics, especially those requiring numerical derivatives, are nontrivial to calculate in 
the SMMC, but in Ref. j25( a method was found to accurately calculate the heat capacity in the temperature range of 
interest. A slight enhancement over the Fermi gas heat capacity was found in even iron nuclei around T ~ 0.8 MeV. 
This enhancement is more pronounced for the neutron-rich iron isotopes and is correlated with the reduction in the 
number of J = neutron pairs with temperature. However, because the calculation was restricted to a finite space 
(fpgs/2 shell), the calculated heat capacity reached a maximum around temperatures of ~ 1.5 MeV and the effect of 
the pairing phase transition was less pronounced. 

Here we use the extended theory to display the heat capacity in a much broader temperature range (i.e., up to 
~ 4 MeV), where deviations from the basic linear dependence become obvious. The effect of the pairing transition 
on the heat capacity is now clearly observed even in the iron isotopes with a smaller number of neutrons in which 
the effect was not noticeable before, e.g., 56 Fe. Fig. [5] shows the extended heat capacity of 56 Fe calculated from 124|) 
and H34f> (solid squares) in comparison with the independent-particle heat capacity calculated in the full space (solid 
line). Contrary to the SMMC heat capacity, the extended heat capacity continued to increase monotonically with 
temperature as does C sp . However, it also exhibits a 'shoulder' at low temperatures that is a clear signature of the 
pairing phase transition and is qualitatively similar to the measured heat capacity in even-even rare-earth nuclei |24j . 



V. HEAT CAPACITY AND THE PAIRING PHASE TRANSITION 



12 



40 



30 



O 20 



10 




1 2 3 

T(MeV) 



FIG. 9: The extended heat capacity of Fe versus temperature, showing a bump around the pairing transition temperature 
(solid squares). The SMMC heat capacity is shown by the open squares. The dashed line is the heat capacity obtained from 
the independent-particle model (bound states plus continuum). 



VI. CONCLUSION 



In this work we have combined the partition function of the independent-particle model in the full space with the 
partition function of the interacting shell model calculated in a small space by the SMMC method. This enables us to 
extend the SMMC calculations of partition functions and level densities up to significantly higher temperatures. The 
results for iron nuclei suggest that the back-shifted Bethe formula i)27|) (with temperature-independent parameters) 
is valid over a wide range of temperatures extending from T ~ 1 MeV to T ~ 4 MeV. The two parameters of the 
formula, the single-particle level density parameter a and the back shift A, agree well with the empirical values. 
A corresponding empirical formula l|19|l that accommodates correlation effects describes well the logarithm of the 
excitation partition function. 

Expressed in the form a = A/K, the theoretical value for 56 Fe is K — 9.5 MeV. This is considerably smaller than the 
Fermi-gas value of K « 15 MeV. It is also smaller than the value K « 10.7 MeV extracted from the high temperature 
slope of the logarithm of the independent-particle partition function, indicating the importance of correlation effects. 
In the literature, effective values of the parameter a are often defined from relations that are valid in the Fermi gas 
limit (and that ignore corrections due to particle- number fluctuations); e.g., E x — aT 2 , S — 2aT and S* 2 = iaE 2 . In 
the presence of shell and correlation effects, these effective values of a usually differ from each other and exhibit a 
temperature dependence even at lower temperatures. Using the empirical relation (|19fl for the partition function, we 
find a constant value of a that is similar to the value extracted from the level density itself. 

The effects of interactions and unbound states were also taken into account in Ref. |l2j| . The authors parameterized 
their results with a temperature-dependent K, and found it to vary significantly, unlike our K . However, their work 
only included part of the interaction, omitting in particular the pairing interaction (which is largely responsible for the 
back shift), and their treatment of the continuum was approximate. The variation they observe at low temperatures 
T ~ 1 — 2 MeV is at least partly due to definitions of K that are based on the Fermi gas limit. 

For temperatures below ~ 1 MeV, the pairing effects are strong and beyond the range of applicability of simple 
parameterizations. This can be seen in the canonical heat capacity if it is displayed over a large temperature range, 
as we did using our extended theory. 
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APPENDIX A 

In this Appendix we derive the empirical back-shifted formula (|19f) for the excitation partition function. The 
partition function is obtained from the level density by a Laplace transform 



Z'(p) = / p{E x )e-P E *dE x . (28) 
Jo 

Substituting the back-shifted Bcthe formula (|27Jl for p{E x ) in Eq. Il28[l . and defining x — E x — A, we have 



We evaluate the integral in l|29|) in the saddle point approximation. The saddle point is given by xq — aT 2 , and we 
find 



APPENDIX B 

In this Appendix we obtain explicit expressions for the canonical energy and heat capacity in the independent- 
particle approximation in terms of logarithmic derivatives of the grand-canonical partition function. 

We introduce the following notation for the logarithmic derivatives of the grand-canonical partition function in the 
independent-particle model 

* ab ~ da"d(J b • [6 J 

Relation © determines a = a({3, N). Considering a a function of (3 and N in the relation (JSJ, we find from dN/d(3 = 

da Z\ i 



d(3 Z 2fi 

We can now use Eqs. {HJ and ill'L'l) to obtain 



(32) 



E sp - -Z ,i + Tf^- o'yz 1 ' 1 ■ ( 33 ) 

■">2,0 -^2,0 

Similarly, the canonical heat capacity in the independent-particle model is given by 

2 ( „ 22^+^2,2 -2^2 1+2^3^^1.1+23.0^1,2 ^Zz$Zl,\Z\ y \ + ZiflZ^ Z\qZ\ x \ 

C W ~ ~P ~ Z 0,2 + ^ ^2 + ^3 • ( d4 ) 

\ >">2,0 Z ^2,0 Z ^2,0 ^2,0 / 

Eqs. (|33|l and (|34|l are used to calculate the canonical thermal energies and heat capacities in the independent- 
particle model in both the full and truncated spaces. 

APPENDIX C 

In this Appendix we derive explicit expressions for the partial derivatives Z a ^, of the grand-canonical partition 
function (see Appendix B) in terms of the single-particle spectrum and scattering phase-shifts. 
Using Eq. we find 



lj 



Snij) b g a +b(a - /3e n ij) + — [ de^-(-e) b g a „ 
7T J de 



b (a - fie) 



(35) 



where we have defined 

^)=ln(l + e-); g a = £§ . 
Note that gi(x) — 1/(1 + e x ) = f, where / is the Fermi-Dirac occupation number for x — a — fie. We also have 



.92 


= /l 


= /(l 


-/) 






93 


= /a 


= /(l 


-/)(1 


-2/) 




94 


= h 


= /(l 


-/)(1 


-6/4 


6/ 2 ) 


95 


= h 


= /(l 


-/)(1 


-14/ 


+ 36/ 2 - 24/ 3 ) 



To avoid the numerical derivative of the phase-shift, we can integrate l|35|) by parts to obtain 

{b f°° 
<iij9a+b{a - fJe nlj ) - - / de S h (e)e b - 1 g a+b (a - fie) 



+ — I de Sij(e)e b g a+b+ i(a - fie) 



for b > and 

ZaS) = X]( 2 -? + \ \9a{a - /3erUj) ~ 9a(a)} + - ( de 5 lj (e)g a+1 (a - fje) 

for b = 0. 

We note that Eqs. I jlljl and (|12|l can be rewritten as 

N - Z 0>1 ; ((AiV) 2 ) = Z ,2 , 
where Zq^ and Zo,2 are given as special cases of (|38p. 
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